Centro-symmetric Hamiltonians foster quantum transport 
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We propose a model for fast and highly efficient quantum transport of excitations, through finite, 
disordered systems. The presented mechanism is statistically robust against configurational changes 
which alter the realization of disorder. We finally discuss the potential relevance of our findings for 
excitation transport in photosynthetic light harvesting complexes. 
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In a variety of fields, ranging from quantum informa- 
tion [l[ to solar cell physics 0, the efficient transport of 
quanta is of paramount importance. Under which general 
conditions fundamental principles of quantum mechanics 
can be exploited to enhance transport in complex systems 
remains a widely open question. Common wisdom sug- 
gests that quantum interference can enhance transport 
across perfectly periodic potentials 0, H| , while it tends 
to suppress transport in disordered systems [H, @]. In 
general, multi-path quantum interference leads to erratic 
fluctuations of transmission probabilities when boundary 
conditions are slightly changed [7- 11 1. These fluctuations 
are often indicative of the strong, non-linear coupling of 
few degrees of freedom, as it abounds in heavy nuclei 
Q , ultra-cold many-particle dynamics [l2j , strongly per- 
turbed Rydberg systems 3;15| billiard geometries for 
photons 1 1611 and electrons 17 1 , strongl y d riven quantum 
systems |l8| . and in large molecules 1|| 2^|. Can we iden- 
tify structural properties of the underlying Hamiltonians 
in those specific cases which exhibit quantum- enhanced 
transport? Are these properties statistically robust, in 
the sense that they are "implementable" by controlling 
only few coarse-grained parameters, without claiming full 
control of the detailed structure of a possibly large, com- 
posite quantum system? It is the present contribution's 
purpose to offer an affirmative answer to these questions, 
on the basis of very general considerations. As a possi- 
ble application, we will discuss the potential role of our 
findings for efficient light harvesting in photosynthetic 
complexes. 

Let us start with the definition of our model, which is 
inspired by the structural elements of the light harvest- 
ing complex FMO of sulfur bacteria 2l|, |22| and by the 
theory of quantum graphs We consider the coher- 
ent quantum transport of one single excitation across a 
disordered network of N sites. Hilbert space is spanned 
by the basis states \i) which represent those states where 
the excitation is fully localized at the network's site i. 
In order to formulate a quantitative, statistical theory, 
we generate different realizations of disorder by sampling 



over N x N Hamiltonians H from the Gaussian Orthogo- 
nal Ensemble [24j, with matrix entries Hij which encode 
the couplings between sites i and j. For each realization, 
input | in) and output |out) are defined as those sites with 
the weakest coupling V = minj^y \Hij\. Our figure of 
merit is the transfer efficiency 



V H = max |(out,0(t))| 2 , |0(O)) 
te[o,T H ) 



(1) 



which quantifies a given random network's performance 
in terms of excitation transport from |in) to |out). Vh 
is gauged against the direct coupling V between |in) to 
|out), in the absence of all intermediate sites, through the 
definition of the associated benchmark time scale = 
tt/2V [13, [25]. 

We have shown earlier, in a study of 3D random graphs 
with dipole-dipole coupling j2(| that a strong and unam- 
biguous correlation exists between high transfer efficien- 
cies and the Hamiltonian's vicinity in matrix space to a 
centro-symmetric structure with respect to |in) andjout), 
which is defined by JH = HJ and |in) = J|out) [27j]. J 



is the exchange matrix, 



Si. 
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However, we 



did also observe that even for high ccntro-symmctries the 
distribution of transfer efficiencies was still rather broad, 
implying that centro-symmetry alone is not sufficient for 
efficient state transfer. This can be easily understood 
by the fact that, while centro-symmetry does distinguish 
| in) and |out) from the intermediate sites, it does not ex- 
clude network conformations which trap the excitation in 
the vicinity of |in), and thus inhibit transmission. There- 
fore, we need to identify an additional structural clement 
which guarantees robustness, in the sense that the trans- 
fer efficiency must not depend strongly on the specific 
conformation of the intermediate sites. Note that such a 
feature is also of obvious relevance for our model to be 
applicable to realistic light harvesting complexes, which 
indeed do undergo structural changes on the macromolec- 
ular scale. 

Intuitively, this structural stability of efficient excita- 
tion transfer from |in) to |out) is guaranteed if both states 
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are coupled through a dominant tunneling doublet in the 
spectrum. The sole role of the intermediate states is then 
to collectively amend the effective tunneling coupling by 
a collective energy shift As. If As has the proper sign, 
this can lead to a dramatic enhancement of the transfer 
efficiency. Note that such collective shifts induced by the 
coupling to random or "chaotic" states have been inves- 
tigated in the context of chaos assisted tunneling (CAT) 
23-3l|. The existence of a dominant doublet will en- 
ter as a key hypothesis into the subsequent analytical 
description of our problem. 

Given the centro-symmetry of H, it can be cast, 
through an orthogonal transformation to the eigenbasis 
of the exchange operator J, into the block diagonal rep- 
resentation [281 



H = 



H+ 
H~ 



(2) 



In this new form, both H + and H~ are again N/2 x N/2 
GOE matrices, i.e. the elements are sampled from a 
Gaussian distribution with zero mean and variance (1 + 
6 h3 )2e/N. 

Since two of the eigenvectors of J have the form 



\±) 



V2 



(|in)±|out)), 



(3) 



we now additionally assume (see above) that |+) and |— ) 
form a dominant doublet, such that they are both close 
to eigenstates |+) and | — ) of H + and H~ , respectively. 
It is then useful to write the Hamiltonian ([2]) as 



H = 



(E + V (V+| 
|V+) Hf ub 

\ 



\ 

E-V (V~\ 
IV") H-J 



(4) 



which makes the definition of rows and columns which 
relate to |+) and |— ) explicit. From the definition of 
|±) it is easy to see that (H)± = E ± V. |V ± ) contains 
the couplings of the dominant doublet states |±) to the 
remainder of the system. 

Since the norm HV^I of the coupling is small, due to 
the assumption of the dominant doublet, 



K±-±)l 



> a 



1 



(5) 



Under this condition, perturbation theory guarantees 
that E ± V in Q are eigenvalues of H , up to some small 
correction term s^. The explicit expression for the trans- 
fer efficiency is then dominated by those terms associated 
with |±), leading to the estimate 



Vh ^ max — e 
te[o.T R ) 4 I 



-it(E+V+s + ) 



j-it(E-V+s~ 



where s = 



E 



E ± V - ef ' 
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and \ip^) and ef are the eigenvectors and eigenvalues 
of H^ ub . From ^ it is clear that the efficiency is large, 
Vh > ct 2 -, if t = tt/\2V + As| is smaller than Tr, and we 
can interpret \2V + As| as an effective tunneling rate. All 
realizations for which we obtain Tji/t > 1 have efficient 
transport which is faster than the direct coupling between 
| in) and |out). 

For fixed E and V the distribution of As was already 
derived within the context of the CAT [3(J |3l|, and it 
turned out that for large N it is a Cauchy (or Lorentzian) 
distribution. In our present problem, E and V are 
stochastic variables, and we should average the resulting 
distribution over E and V. Considering that the integra- 
tions over E and V arc dominated by the mean values of 
their distributions, which are given by V w -y/^^/A^ 3 / 2 
and E 2 = 2£ 2 /N (since E = 0, we use E 2 ), a lengthy 
but straightforward calculation shows that the probabil- 
ity distribution of Tr(2V + As) /it = Tn/t is given by 



with sq 



f 



so 



+ 



so 



ir \s 2 + (I + x) 2 s 2 + (l-a;) 2 
jVj|2ArV2(2-i/2A0 



2tt£ 2 
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Here, ||V|| 2 denotes the expectation value of UV^I 2 for 
all realizations where the dominant doublet assumption 
holds. 

Under the prerogative of a dominant doublet, the dis- 
tribution (0 depends on only two parameters: £ charac- 
terizes the spectral density of the eigenstates of H* ub and 
H~ ub , while || V|j 2 measures the average coupling strength 
of the dominant doublet to these states. It therefore can- 
not be emphasized enough that, within the picture here 
elaborated, the transport properties of the problem do 
not depend on the specificities of the intermediate elec- 
tronic states of the network - i.e. on its specific structural 
conformation. 

In order to validate our theoretical model by numer- 
ical simulations, we generate many GOE Hamiltonians 
with the additional constraint of centro-symmetry with 
respect to |in) and |out). For each of these Hamiltoni- 
ans the existence of a dominant doublet is assessed by 
inspection of its eigenvectors and by verifying that, for 
some |±), the condition |(±,±)| > a holds. Then, the 
quantity Vh is obtained by numerically solving the quan- 
tum dynamics of the Hamiltonian, and t is defined as the 
earliest point in time for which |(out, <f>(t)) | 2 = Vh- 

Figure [T] compares numerical results and analytical 
prediction (J7J. It must be noticed that there are no 
free parameters involved in this plot: The values £ = 2, 
a > 0.95 and ||V|| 2 = 0.323192 are either given a pri- 
ori or directly extracted from the statistical characteri- 
zation of the numerically generated sample Hamiltonians. 
Clearly, the majority of realizations have t smaller than 
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FIG. 1: Histogram of the numerically simulated inverse trans- 
fer time Ta/t (dashed line), compared to the theoretical pre- 
diction (0 (full line), for density of states £ = 2, N = 10 
network sit es, an d a dominant doublet condition a — 0.95 
(see (©). ||V|| 2 = 0.323192 is directly extracted from the 
numerical sampling over different random Hamiltonians, and 
the evident agreement between simulation and theory results 
with no adjustable parameters. Simulations were performed 
in the time window [0, 1.7Tr]; therefore the plot only starts 
at Tn/t — (1.7) . The distribution's heavy tail (as high- 
lighted in the inset, on a log scale) guarantees the abundance 
of very short transfer times t <C Tr as compared to the direct 
coupling time between |in) and |out). 



the time scale set by the direct coupling V between |in) 
and | out). The fat algebraic tail of the Cauchy distribu- 
tion for t << Tji guarantees that realizations with very 
fast transport, with respect to Tr, are abundant in the 
sense that they are not exponentially unlikely. 

We still need to verify that the such defined Hamilto- 
nians provide efficient transport. For a fixed a = 0.95, 
our model predicts that the efficiency should be larger 
than Vh ^ c? « 0.9. This is indeed observed in the sim- 
ulations. Figure [2] shows the probability distribution of 
the efficiencies, which is sharply peaked above Vh > 0.9. 
Comparison, in the same figure, with the probability dis- 
tributions of the efficiencies of centro-symmetric GOE 
matrices without the doublet constraint, and of general 
GOE matrices, respectively, shows that in both cases 
the average efficiency is significantly lower than for these 
ccntro-symmctric Hamiltonians which exhibit a domi- 
nant doublet. 

A remarkable asset of the transport mechanism here 
presented is its robustness under different realizations 
of disorder, which, in the context of networks, refers to 
different configurations of the intermediate sites (repre- 
sented by H^ h in (|4|)). In the light of the recent debate 
on the potential role of quantum coherence in photosyn- 
thetic harvesting of the sunlight's energy, one may won- 
der whether this mechanism is already implemented by 
Nature. We now briefly scrutinize such hypothesis, on 



FIG. 2: Histograms of the transfer efficiencies Vh for three 
different Hamiltonian matrix ensembles, with fixed density of 
states £ = 2 and number iV = 10 of network sites. Clearly, as 
we progressively constrain the ensemble from GOE to centro- 
symmetric GOE and, finally, to centro-symmetric with dom- 
inant doublet condition a = 0.95, the average efficiencies are 
dramatically enhanced. In perfect accord with our analytical 
model, the efficiency distribution is strongly peaked at val- 
ues Vh ^ o? ~ 0.9 (indicated by the arrow), when we only 
sample over random Hamiltonians with dominant doublet. 



the basis of published data on the widely studied FMO 
light harvesting complex. 

To start with, the recently discovered eighth site of 
the FMO molecular network is only weakly coupled to 
the remainder of the complex plj . This is indeed con- 
sistent with the above intuition of a dominant dou- 
blet. Furthermore, the centro-symmctry of the pub- 
lished FMO Hamiltonian H 21] was studied in pfjj ]. 
and it was found that, while Hq is itself neither effi- 



cient nor centro-symmetric, weak perturbations of the 
matrix elements of Hq allow for high efficiencies [HJ , to- 
gether with enhanced centro-symmetries (26|. Consis- 
tently, verification of the dominant doublet condition ([5]) 
for Ho [33| results in max^l (±, r]i) | 2 ~ 0.5, upon maxi- 
mization over all eigenvectors 1 77^) . In contrast, the op- 
timized Hamiltonian inferred in [26| indeed does exhibit 
a dominant doublet, with max^ |(+,»7i)r = 0.877361 and 
max, 

\(-,m)\ = 0.901908, respectively. 
In conclusion, we described a general mechanism that 
gives rise to fast and efficient quantum transport in fi- 
nite, disordered systems. The mechanism rests on two 
crucial ingredients: The first is the centro-symmetry of 
the underlying Hamiltonian, which guarantees the exis- 
tence of a block diagonal representation. The second is 
the existence of a dominant doublet, that guarantees the 
statistical robustness of the transport properties, under 
the coupling to random/chaotic states which assist the 
transport. Note, in particular, that the structure of (3]) 
is maintained under averages over different realizations 
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of centro-symmctric matrices with dominant doublets 
(what provides an explanation for the on a first glance 
surprising numerical observation in pril ]. that also aver- 
ages over efficient Hamiltonians remain efficient). Both 
these features appear consistent with published effective 
Hamiltonians of the FMO light harvesting complex and 
weak perturbations thereof. The statistics of the transfer 
efficiencies and times as shown in Figs. 1 and 2 only de- 
pend on the intermediate network sites' densit y of states 
£, and on the average coupling strength ||V|| 2 of the in- 
and output-sites to the network. On the one hand, this 
may explain why coherence effects can survive even in 
relatively floppy macromolecular structures, if only these 
average quantities can be stabilized by the functional 
units' architecture. On the other hand, if such stabi- 
lization is possible, this immediately implies a handle of 
control, by proper tuning of £, ||V|| 2 and, possibly, TV 
21| . according to the physiological needs. 
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